SIGNAL-TO-NOISE RATIO TUNED ADAPTIVE OPTICS CONTROL SYSTEM 



CROSS REFERENCE TO RELATED APPLICATIONS 
None. 

5 STATEMENT REGARDDsfG FEDERALLY SPONSORED RESEARCH OR 
DEVELOPMENT 
None. 

BACKGROUND OF INVENTION 

The present invention relates to an improved adaptive optics control system and in 

10 particular, to a signal-to-noise ratio-tuned wavefront corrector for use therein. 

Optical wavefronts are subject to distortion when passed through certain mediums. 
These distortions can degrade the quality of images of an object being observed through an 
image-forming device such as a camera or telescope. Such distortions and degradations can be 
especially severe and problematic in turbulent mediums such as the atmosphere. Some of the 

15 diffractive effects of such turbulence include image blurring caused by the turbulence 
aberrations, scintillation caused by propagation through strong turbulence, and speckle caused by 
coherent scattering from diffuse objects. 

Various adaptive optic control systems, mechanisms and methods have been developed to 
try to correct for the blurring that results from such distorted wavefronts. One such control 

20 system for use with telescopes, often referred to as a "tip-tilt" correction, involves tilting the 
secondary mirror of the telescope several times a second to reduce or eliminate the dancing 
motion of the image. This method, however, only provides a small improvement in the 
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sharpness of the image in applications involving large telescopes. 

Various other adaptive optics control systems have been developed to try to better 
compensate for such distortions. Examples of various such systems are disclosed in "Adaptive 
Optics for Astronomical Telescopes," by John W. Hardy, 1998 (pages 31-33, and 55-69), which 

5 is incorporated herein by reference. Such control systems typically include a wavefront slope 
sensor for measuring the phase differences or phase slopes between points of a wavefront, a 
wavefront reconstructor for estimating the wavefront phase from the phase differences, a control 
system for reducing the effects of noise, and a wavefront corrector for correcting the wavefront 
based thereon. The wavefront sensor is usually in the form of a Hartmann wavefront sensor as 

10 shown in Fig. 2. Hartmann wavefront sensors use an array of lenslets or a mask pierced with 
holes for dividing the distorted wavefront into an array of subapertures. Each of the beams in the 
subapertures is focused onto one or more detectors disposed behind the holes or lenslets. When 
the distorted wavefront passes through the holes or lenslets, it forms an array of spots on the 
detectors, which are indicative of the local wavefront slope or tilt, if any, in the corresponding 

15 subaperture. Typically, the wavefront slope sensor includes an analog-to-digital converter and 
one or more processors to compute the wavefront slopes. 

The wavefront corrector is usually in the form of a deformable mirror which compensates 
for such distortions. Deformable mirrors typically comprise a face plate to the back of which a 
plurality of actuators are secured. The actuators expand or contract in length upon application of 

20 a voltage or a magnetic field in accordance with the electrical conmiands generated by the 
wavefront reconstructor, thereby pushing or pulling on the faceplate and causing the mirror to 
change its shape to make the appropriate corrections to the distorted waveform. The actuators 
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are typically arranged in a square or hexagonal array defining a plurality of zones, and are 
capable of displacing the faceplate locally within each zone by a few micrometers up or down. 

Current wavefront slope sensors, however, are subject to various measurement errors that 
degrade wavefront correction. First, there is noise in each measurement, including read-out 

5 noise and dark-current noise in the detector and shot noise in the received light. Second, 
diffractive cross-talk can be produced among adjacent subapertures. Third, the subapertures and 
detectors may be misaligned. Fourth, partially filled lenslets can cause stretching of the spot size 
on the detectors. Fifth, the intensity across each lenslet may not be uniform, which can also 
change the spot sizes and shapes. In addition, Hartmann wavefront sensors are unable to provide 

10 accurate wavefront phase slope measurements in some subapertures during conditions of severe 
turbulence due to large variations in the intensity and phase within a subaperture. One approach 
to minimize such effects is to combine a Hartmann sensor with another device, such as a unit 
shear lateral shearing interferometer wavefront sensor as is disclosed in U.S. Patent Number 
4,518,854, which is incorporated herein by reference. 

15 Further, conventional control systems assume that all of the subaperture-based 

measurements are of equal quality and that all of the computed signals for the wavefront 
corrector actuators are of equal quality. However, the measurements in the edge subapertures 
(i.e., around the outer aperture, around interior obscurations, or in the vicinity of struts) are 
typically only partially filled, and so have lower signal-to-noise ratios and are of a lesser quality 

20 than the measurements in the interior subapertures. The wavefront correction signals for edge 
actuators are also typically of a lesser quality than those corresponding to non-edge actuators, 
due in part to the fact that the quantity of data in the vicinity edge actuators is limited and due to 
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the fact that the estimate of the phase for each actuator is dominated by the data in the vicinity of 
the corresponding subaperture. Therefore, the correction signals for such edge actuators can be 
extremely noisy, making the system unstable. As a result, the overall bandwidth of such control 
systems is often reduced to the level of the least stable actuator estimate. Such a limitation on 

5 the bandwidth, however, can significantly degrade the performance of the system. 

As a result, there is a need for an improved, adaptive optics control system which 
overcomes these problems. 
BRIEF SXJMMARY OF THE INVENTION 

An adaptive optics control system for distorted optical wavefronts having a wavefront 

10 corrector, a wavefront slope sensor, a wavefront reconstructor and a wavefront controller is 
disclosed. The wavefront corrector has a surface controlled by a plurality of actuators. The 
wavefront slope sensor has subaperture separation means for defining a plurality of subapertures 
through which the distorted wavefront can pass, the wavefront slope sensor being adapted to 
measure the wavefront slope through each subaperture and generate a wavefront slope sensor 

15 output signal for each subaperture indicative of the distortion of the wavefront. The wavefront 
reconstructor is adapted to receive the wavefront slope sensor output signals and calculate a 
phase estimate for each of the actuators of the wavefront corrector, the phase estimate having a 
signal-to-noise ratio. The wavefront reconstructor is further adapted to generate a plurality of 
correction signals, each having a bandwidth, based on said phase estimates. The wavefront 

20 controller is adapted to selectively adjust the bandwidth of each correction signal based on the 
signal-to-noise ratio of the estimated phase calculated by the wavefront. 

In one embodiment, the wavefront controller is adapted to reduce the bandwidth of the 
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correction signals for actuators having corresponding phase estimates with a low signal-to-noise 
ratio, and increase the bandwidth of the correction signals for actuators having corresponding 
phase estimates with a high signal-to-noise ratio. In another embodiment, the bandwidth of each 
correction signal is proportional to the signal-to-noise ratio of the corresponding estimated phase. 

5 The wavefront corrector may comprise, without limitation, a deformable mirror or a spatial light 
modulator. The subaperture separation means may comprise, without limitation, an array of 
lenslets and a corresponding detector array having a plurality of pixels disposed behind each 
lenslet, each lenslet defining a subaperture, a lateral shearing interferometer, or both. The 
wavefront slope sensor preferably comprises means for individually calibrating each wavefront 

10 slope measured by the wavefront slope sensor. 

A method of optical wavefront distortion correction using a wavefront corrector having a 
surface controlled by a plurality of actuators is also disclosed. The method comprises measuring 
the distortion of the wavefront with a wavefront slope sensor having subaperture separation 
means for defining a plurality of subapertures through which the distorted wavefront can pass, 

15 the wavefront slope sensor being adapted to measure the wavefront slope through each 
subaperture and generate a wavefront sensor output signal for each subaperture indicative of the 
distortion of the wavefront, calculating a phase estimate for each of the actuators of the 
wavefront corrector based on the wavefront sensor output signals, each phase estimate having a 
signal-to-noise ratio, generating a plurality of correction signals to be applied to each actuator 

20 based on the phase estimates, each correction signal having a bandwidth, and selectively 
adjusting the bandwidth of each correction signal based on the signal-to-noise ratio of the 
corresponding estimated phase. The method may further comprise reducing the bandwidth of 
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correction signals of actuators having corresponding phase estimates with a low signal-to-noise 
ratio, and increasing the bandwidth of correction signals of actuators having corresponding phase 
estimates with a high signal-to-noise ratio. The method may further comprise selectively 
adjusting the bandwidth of each correction signal to be proportional to the signal-to-noise ratio of 
5 the corresponding estimated phase. The method may also further comprise individually 
calibrating each wavefront slope measured by the wavefront slope sensor. The step of 
calculating a phase estimate may comprise representing the wavefront corrector and wavefront 
slope sensor measurements as a plurality of linear equations in a matrix format. 
BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWINGS 
10 Fig. 1 is a block diagram representing an adaptive optics control system in accordance 

with one embodiment of the present invention. 

Fig. 2 is a diagram of one embodiment of the wavefront sensor of Fig. 1. 

Fig. 3 A is a diagram of the detector of Figure 2. 

Fig. 3B is a diagram of a detector zone of the detector of Fig. 3 A. 
15 Fig. 4 is a cross-sectional view of one embodiment of the wavefront corrector of Fig. 1. 

Fig. 5 is a diagram of the signal processing of the adaptive optics control system of Fig. 

1. 

DETAILED DESCRIPTION OF THE INVENTION 

An adaptive optics control system 10 for correcting distortions in a wavefront 12 in 
20 accordance with one embodiment of the present invention is shown in Fig. 1. While the system 
will be discussed in connection with a telescope 14 for space observations, it can be appreciated 
that it can be used with any device and in connection with any medium involving optical 
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wavefronts which are subject to distortion. The system 10 is a closed loop system consisting 
generally of a wavefront slope sensor 16, a wavefront reconstructor 26, a wavefront control 
system 27 and a wavefront corrector 28. 

While the wavefront corrector will be discussed in connection with a deformable mirror, 

5 it can be appreciated that any spatial light modulator may be used and the term deformable 
mirror as used herein applies to both. As shown in Fig. 4, the deformable mirror 28 comprises a 
thin reflective surface 30 controlled by a plurality of actuators 32 secured to the back thereof. In 
a preferred embodiment, the actuators 32 are piezo-electric actuators, but it can be appreciated 
that any type of actuator may be used. The size of a deformable mirror may vary depending on 

10 the application, ranging in size from about five centimeters in diameter for a small mirror with 
21 actuators, to up to about forty centimeters in diameter for a large mirror with 2,000 actuators. 
In addition, while not shown, a steering mirror may be positioned in front of the wavefront 
corrector in order to optimize the correction range of the deformable mirror. 

Referring back to Fig. 1, the wavefront slope sensor 16 receives the distorted wavefront 

15 12 and measures the departure of the wavefront from a plane wave. One well known 
embodiment of a wavefront sensor known as a Hartmann wavefront sensor is shown in Fig. 2 
and disclosed in U.S. Patent No. 4,518,854. However, other sensors may be used, such as a 
lateral shearing interferometer as disclosed in U.S. Patent No. 6,163,381, and Shack, Roland V. 
and B.C. Piatt, "Production and use of a Lenticular Hartmann Screen," J. Opt. Soc. Am. (Journal 

20 of the Optical Society of America) 1971, pages 67, 81-86, both of which are incorporated herein 
by reference. The wavefront slope sensor 16 uses an array of lenslets 18 for dividing the 
distorted wavefront 12 into a plurality of subapertures. It can be appreciated, however, that any 
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means for subaperture separation may be used. Wavefront slope sensor 16 includes a 
corresponding detector array 20 having a plurality of pixels (not shown) disposed behind each 
lenslet 18 of the wavefront slope sensor 16 and aligned therewith. Typically, the detector array 
is a focal plane array, such as a charge coupled device (CCD) camera. However, any type of 

5 detector array may be used. 

In a preferred embodiment as shown in Figure 3 A, the array of lenslets 18 and detector 
array 20 are each in the form of a two-dimensional array or matrix made up of a plurality of 
subaperture zones (not shown) and detector zones 24, respectively. The size of each array may 
vary depending on the particular application. The number of lenslets 18, however, is preferably 

10 selected to correspond to the number of actuators 32 of the wavefront corrector 28. As shown in 
Fig. 3B, when the distorted wavefront 12 passes through the lenslets 18, it produces a 
corresponding array of spots 22 within each detector zone 24. The location of the centroid of 
each spot 22 in the detector array 20 is indicative of the local wavefront slope or tilt, if any, in 
the corresponding subaperture. In a preferred embodiment, separate wavefront sensors are used 

15 to sense tilts in two orthogonal axes; namely the X and Y axes. However, a single wavefront 
sensor may sense tilts in both orthogonal axes simultaneously. 

In a preferred embodiment, the slope of each subaperture is computed measuring the 
centroid of the light on the pixels of each subaperture and multiplying this centroid by a constant 
to convert it to a slope. The centroid may be computed in a number of ways including, without 

20 hmitation, using all of the subapertures pixels, or a predetermined subset, or using just the pixels 
with values above some threshold. The slope conversion constant depends on the focal length of 
the lenslet, the wavelength of the light, the size of the pixels, and on the units needed for the 
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wavefront reconstructor 26 (e.g., radians of tilt or delta phase per subaperture). The slope 
conversion factors are calibrated for each subaperture by applying known tilts to the wavefront. 
The ratio of actual tilt to the measured slope minus the bias is the slope conversion factor. The 
bias is represented as the slope measured when an undistorted wavefront is detected. The biases 
5 can be computed in several ways. In one embodiment, a local laser with a clean (flat) wavefront 
is put into the system. This local signal is averaged over many measurements (frames) to 
remove the effect of noise. A few hundred frames is generally sufficient. Since the input signal 
is error-free, the averaged measured signal is the bias. In another embodiment, a reasonably 
bright star is imaged and the measurements are averaged over many seconds of data to average 
10 out not only the noise, but also the atmospheric turbulence. With such a configuration, the bias 
and slope-conversion factors can be individually calibrated for each subaperture of the wavefront 
slope sensor 16. The output slope computed by the wavefront slope sensor 16 is then the "raw" 
slope minus the bias. 

The output of the wavefront slope sensor 16 consists of a plurality of output signals 
15 corresponding to an array of zonal measurements from each detector zone 24 representing the 
local wavefront slopes of each subaperture in two orthogonal axes. Referring back to Fig. 1, the 
wavefront sensor output signals are transmitted to the wavefront reconstructor 26 which then 
generates a corresponding plurality of wavefront phase estimates based thereon. The wavefront 
reconstructor 26 includes an analog-to-digital converter (not shown) for converting the analog 
20 output signals of the wavefront slope sensor 16 to digital output signals, and one or more 
processors (also not shown) for processing the digital output signals to generate the correction 
signals needed to correct the distorted wavefront. It is understood by one skilled in the art that 
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many processors are suitable for this purpose. 

The wavefront reconstractor 26 translates the subaperture tilt measurements represented 
by each digital output signal into a plurality of phase estimates for each actuator of the 
deformable mirror 28. Specifically, given the measured slopes sensed by the wavefront slope 

5 sensor 16, each having some measurement error, the wavefront reconstructor 26 determines the 
best estimate of the wavefront phase at the actuator locations of the wavefront corrector 28 that 
gave rise to the slopes measured by the wavefront slope sensor 16. There are multiple 
approaches to wavefront phase reconstruction, any one of which may be used without affecting 
the innovation of the present invention. One approach constructs multiple paths between 

10 actuators and averages them. Another approach uses well known statistical methods. The most 
conmion method is to treat the wavefront phase reconstruction process as a system of linear 
equations represented in a matrix format. The mathematics associated with this approach is well 
known to those skilled in the art of numerical methods for linear algebra. One method is 
described in William H. Press, et al. Numerical Recipes in C: The Art of Scientific Computing , 

15 Cambridge University, 1992, which is incorporated herein by reference, and routines for this 
method are available in software such as The Math Works, Inc.'s MATLAB product. 

One embodiment of solving the wavefront phase reconstruction problem is as a system of 
linear equations, defined by the following parameters: 
M = the number of actuators; 

20 N = the number of subapertures (typically, N = 2M); 
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L = the vector of wavefront slopes in the x-direction measured by the wavefront 

slope sensor; 



= the vector of wavefront slopes in the y-direction measured by the wavefront 
slope sensor; 



s = 



= the combined vector of x and y wavefront slopes; and 



= the vector of the heights of the actuators 32 of wavefront corrector 28. 
If an actuator 32 of the wavefront corrector 28 is pushed, then the response to that 
actuator is measured by the wavefront slope sensor 16 and represented in the wavefront slopes. 
A response matrix, /f , of N rows and M columns representing the combination of all the actuators 
32 being pushed on in sequence is defined as follows: 

The expected response to an arbitrary vector, a, of actuator positions is a linear combination of 
the responses to individual actuators. This can be written in matrix notation as: 
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and assumes that the slope response to an individual actuator is hnear with the magnitude of the 
push, and that the response to multiple pushes is the sum of the individual pushes. Since this 
assumption is not perfectly true, it can be improved upon, if required, through experiment by 
pushing on individual actuators and groups of actuators and measuring the surface of the 
wavefront corrector with a commercially available interferometer. 

The solution to the wavefront reconstruction problem (i.e., the estimation of the 
wavefront phase at the actuator locations of the wavefront corrector 28) consists of finding an 
inverse of the actuator response equation set forth above, (equation (1)). A "least-squares" 
approach as is well known to those of ordinary skill in the art is used to find the solution, 2 that 
minimizes the sum of the squares of the differences between the measurements expected from 
the solution and the actual measurements. Namely, p minimizes as set forth below: 



Since R is not a square matrix, p is substituted for a in equation (1) and both sides of the 
equation are multiplied by J?^, the transpose of /?, so that (R^R) is square as set forth below: 



Typically, there is at least one actuator mode or pattern that is unobservable so as to 
create a singular (R^R) matrix. This can occur when all the actuators are poked by equal 
amounts simultaneously, such that the wavefront is still flat and the measured slopes are all zero. 




where 



S = Rp 




(2) 
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The preferred method to solving equations with singular matrices, such as equation (2), is via 
singular-value-decomposition (SVD), as described in section 2.6 of Numerical Recipes in C = 
The Art of Scientific Computing referenced above and which is also available as a package in 
The MathWorks' Inc.'s MATLAB product. SVD decomposes the matrix R^R into other 
5 matrices, which are easier to invert. A summary of this method is as follows: 

UDlf = SVD(R^R) 
where 

1/ is a unitary matrix, so its inverse, IT^ is its transpose, l/; and 

D is a diagonal matrix, with each diagonal element corresponding to a particular mode. 
10 If a mode is singular, its diagonal element will be zero (or nearly zero, due to round-off errors in 
the computer). In such a case, a "pseudo inverse" of D, Z)^, is used whose diagonal elements are 
the reciprocals of the non-singular diagonal elements of D, with zeros for the singular values. 
This causes the singular modes to be ignored. 

Thus, the reconstructor matrix, Wy is computed as follows: 

and the solution to equation (2) is then computed as follows: 
p = WS 

- - (5) 

The negative of g then represents an estimate of the wavefront phase. After the estimate 
of the wavefront phase has been computed, a wavefront controller or control system 27 is used to 
20 filter out the noise. Since the wavefront being measured is usually changing at a slower rate than 
the wavefront slope sensor measurements are being taken, the time sequence of phase estimates 
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can be processed to reduce the influence of the noise on the correction signals. 

The wavefront control system 27 consists of one or more processors, software to 
implement the control algorithm described below, and a digital to analog converter to generate 
the voltages to be applied to the wavefront corrector 28. The wavefront control system 27 is 
5 partially characterized by its bandwidth frequency and by its stability. Wavefront phase 
disturbances of temporal frequencies below the bandwidth are approximately corrected for, while 
phase disturbances of temporal frequencies above the bandwidth frequency are approximately 
ignored. Likewise, the portion of the noise below the bandwidth frequency is passed 
approximately along with the phase correction, and the portion of the noise above the bandwidth 

10 frequency is approximately filtered out. The control system performance is limited by the 
signal-to-noise ratio of the wavefront slope sensor measurements, and by the time it takes the one 
or more processors of the wavefront reconstructor 26 to compute the slopes and to reconstruct 
the wavefront. These factors can lead to instability, so they limit the bandwidth and the 
particular design and performance of the control system. 

15 In the usual implementation of control systems in adaptive optics, the same control 

system equations are applied to each actuator of the wavefront corrector. For many adaptive 
optics applications, however, a low signal contributes to noisy estimates of the phase of the 
wavefront at the actuators. The lowest signal-to-noise ratio actuators are typically edge actuators 
(i.e., those actuators around the outer aperture of the telescope, around interior obscurations, or 

20 in the vicinity of struts) because some of the lenslets/subapertures imaged near these actuators 
are partially obscured, and/or because there are less valid subapertures imaged near such 
actuators. The stability of the control system for these noisy phase estimates then limits the 
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stability of all the phase estimates. This in turn limits the bandwidth and the perfonnance of the 
overall adaptive optics control system. 

In accordance with the present invention, the wavefront control system 27 selectively 
applies a distinct bandwidth to each actuator or correction signal based on the signal-to-noise 
5 ratio of the corresponding estimated phase computed by the wavefront reconstructor 26. In 
particular, the bandwidth is lowered for low signal-to-noise ratio actuators, such as the edge 
actuators, and raised for high signal-to-noise actuators. There may be "typical" signal-to-noise 
ratios between the low and high signal-to-noise ratios in which the bandwidth is neither lowered 
or raised. In one embodiment, low signal-to-noise ratios are less than 5, typical signal-to-noise 

10 ratios range between 5 and 10, and high signal-to-noise ratios are greater than 10. It can be 
appreciated, however, that any signal-to-noise ratio values may be used to classify low, typical 
and high signal-to-noise ratios, and that a different classification scheme other than low, high and 
typical may be used. Also, rather than classifying the signal-to-noise ratios as low, typical, and 
high, the bandwidth of each actuator correction signal may be set proportional to the signal-to- 

15 noise ratio of its estimate, as described below. This allows for an increase in the average 
bandwidth and more filtering of the noisiest estimates, while maintaining or improving stability. 
The performance improvement can be significant, since typically there are only a minority of 
low signal-to-noise ratio actuators. 

The signal-to-noise ratio of each actuator's phase estimate can be computed as follows. 

20 From equation 5 above, the estimate of the phase at each actuator can be written as follows: 

M 

=l,^kiSi (6) 
where Wki is the element in the k^*^ row and i^^ column of the reconstructor matrix, W and Si is i*^ 
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measured slope. Typically, this sum is dominated by the slopes measured in the subapertures 
neighboring the actuator. 

It is well-known in statistics that the variance of a sum is the sum of the variancesand that 
the variance of a coefficient times a random variable is the square of the coefficient times the 
variance of the random variable. The variance of each slope measurement includes shot noise 
(its variance equals the signal) and read noise as set forth below: 



where pi = the signal of the i^^ slope measurement (in photo-electrons) and is the variance of 

the read and dark-current noise, which is usually the same for all subapertures. 

Thus, the signal-to-noise ratio (SNR) for each actuator, k, can be computed as follows: 



The control system for computing the actuator conmiands can be written as a 
combination of the past actuator commands and the current and past reconstructed phase 
estimates, as given by this equation: 



where bk,t = the conmiand for actuator k at time t (counted in iterations); 
at = the phase estimate from the reconstructor for actuator k at time t (counted in 
iterations); and 

Ckj and gk,i are control system coefficients. 

If the control system is a simple integrator that uses only the current reconstmcted phase 




(7) 




(8) 



(9) 
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estimates, then the equation may be written as follows: 

where gk is the control system "gain" for actuator bk. 

The SNR can be used to adjust the gk,i coefficients. The SNR optimized gains are then 
5 calculated as follows: 
SNR, 

where gi is the nominal gain of the control system (not adjusted for SNR) and SNRq is the 
nominal SNR of the wavefront sensor. 

While the invention has been illustrated and described in connection with several 

10 preferred embodiments, it will be appreciated that further changes can be made therein without 
departing from the spirit and scope of the invention. For example, while the invention has been 
discussed with reference to a telescope and the atmosphere, it applies to any device and in 
connection with any medium involving optical wavefronts which are subject to distortion, 
including without limitation any imaging or laser systems. Accordingly, it is intended that the 

15 scope of the invention not be limited by the disclosure of the preferred embodiments, but instead 
be determined entirely by reference to the claims that follow. 
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